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A DETERMINANTAL APPROACH TO IRRATIONALITY 


WADIM ZUDILIN 

Abstract. It is a classical fact that the irrationality of a number ^ G R follows 
from the existence of a sequence Pn/q_n with integral and such that Qn^—pn ^ 
0 for all n and — Pn 0 as n —>■ oo. In this note we give an extension of 
this criterion in the case when the sequence possesses an additional structure; in 
particular, the requirement qn^ — —>■ 0 is weakened. Some applications are 
given including a new proof of the irrationality of tt. Finally, we discuss analytical 
obstructions to extend the new irrationality criterion further and speculate about 
some mathematical constants whose irrationality is still to be established. 


1 . Irrationality criteria 

Let ^ be a real number we wish to prove the irrationality of. Assume there 
exists a sequence Pn/Qn with integral pn and qn such that ^ ^ Pn/qn for all n and 
Qn^—Pn —>■ 0 as n —>■ cxo. If ^ were rational, ^ = p/q say, then q{qn^ — Pn) = QnP—PnQ 
is a nonzero integer for all n, hence its absolute value is at least 1 . On the other 
hand, \q{qn^ ~ Pn)\ is less than 1 for n sufficiently large; a contradiction. 

Let us put the above classical irrationality criterion in a standard context of real 
numbers ^ that happen to be periods. Suppose we have a sequence of rational 
approximations 

rn = an^-bn^Q^ + Q ( 1 ) 

such that 

(a) 0 < < Oig" for some Oi, e > 0 and all n = 1, 2 ,... ; 

(b) SnCin, bnbn ^ ^ for some sequence of positive integers 5 n] and 

(c) 5 n < C'2A" for some O2, A > 0 and all n = 1, 2 ,... . 

Proposition 1 . Under hypotheses (a)-(c), if eA < 1 then ^ is irrational. 

The principal goal of this note is to demonstrate that under some further (natural) 
assumptions on the approximants ([T]) we can replace the condition eA < 1 by a 
slightly different one. Namely, assume additionally that 

(d) = / z{x)'^u{x) for some domain 7 C non-constant continuous func- 

J'y 

tion z{x) > 0 on 7 and measure (positive differential form) uj{x)-, and 
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(e) dn divides 5n+i for all n = 1, 2 ,... . 

Proposition 2. Under hypotheses (a)-(e), if < 1 then f is irrational. 

Our proof, in which the classical Vandermonde determinants as well as Hankel 
determinants of the sequence ([I]) show up, is given in Section|21 Note that eA^/^/4 < 
(eA)3/2 > 1/16, thus Proposition [2] has potentials to produce irrationality results 

when Proposition [1] is not applicable. For example, if A e = 2.7182... and 
0.34 < e: < 0.89 then Proposition |2] implies the irrationality of while Proposition [T] 
does not. In Sections E] and 0] we give such applications of the new irrationality 
criterion to log 3 and vr. Though these numbers are known to be irrational, our 
proofs based on the argument of Proposition |2] are new. 

The example .^ = 1, a„ = 1, &„ = (2”'—1)/2"', so that On^—bn = 1/2”, demonstrates 
the importance of (d): in this case we can take e = 1/2 (or slightly bigger) and 
5n = 2”, hence A = 2. The condition eA^/^/4 = \/2/4 < 1 is clearly satished but 
without any implication, as there is no way to write the linear forms — bn in 
the form assumed in (d) for a non-constant function z{x). (In fact, all the Hankel 
determinants that appear in the proof of Proposition [2] below vanish in this case.) 
On the other hand, the example z{x) = x, u{x) = dx and 7 = [0,1] (this corresponds 
to the choice ^ = 1 or 0), so that e = 1, 

6n = lcm(l, 2 ,... ,n -t-1) 

and A is anything slightly larger than e, shows that the condition eA^/^/4 < 1 
cannot be relaxed “too much.” In Section 0] we discuss the optimality of the latter 
constraint as well as comment on the orthogonality induced by the data from con¬ 
dition (d). Finally, in Section 0] we speculate about some mathematical constants 
whose irrationality is still to be established. 


2. Proof of the irrationality criterion 
Proof of Proposition [21 Conditions (a) and (d) imply that 


0 < supz(a;) < £. 

xG'y 


Consider the polynomial forms 
1 


Rr, = 


n\ 


[■ ■ ■ [ J] (zix^) - z{x^)fu^ e Qr + • • • + Qe + Q. (2) 


l<j<£<n 


Note that Rn > 0, because the integrand is nonnegative. Furthermore, the form 
has a nice Hankel determinant evaluation due to Heine [5] , 

Rn = det (rj+e), 

which together with hypotheses (b) and (e) imply that 


( 3 ) 


dn-ldn • ■ ■ <^2n-2-Rn ^ + ' ' ' + + Z. 
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In addition, we have 

( sup {zj — Zif j e/4 as n ^ oo (4) 

l<7<£<n ^ 

(the Fekete-Chebyshev constant of the interval [0,e]), so that 

Rn < (e/4)”'+°(’"') as n ^ oo. 

If ^ were rational, ^ = p/q, then 11^=0 ^n-i+j ■ q^Rn would be a positive integer 
for any n. On the other hand, 

• • • 52n-2q^Rn < (C'2g)"A(3/2)-^ (e/4)d+°(i))^0 as n ^ oo, 

a contradiction. □ 

In the argument above we could have used the inclusion 11^=0 ^n-i+j ■ q^ • nlR^ G 
Z instead (with the factorial factor) as an immediate consequence of the integral 
representation (j2]), no reference to ([3]) is necessary. But because Heine’s theorem 
about the determinant representation of Rn is itself quite elementary and nice, we 
reproduce its proof here for completeness. 


Proof of identity ([3]) . Using the dehnition of the sequence write 


det (r,_|_£) = det ( / zix^V^^ujixR 

o<j,e<n ^ ^ o<jy<n' I ^ ^ ^ 


det (z(x^y+^) oj(x^) oj(xy • • • u;(x^-y 
0<j,e<n^ 


det u 


0<j,i<n 


^1^2 • • • 4-1 n 

0<j<e<n-l 


Zj)UJ 


where zj = z^x^ and the evaluation 


0<j<i<n-l 


(5) 


is applied. It remains to notice that the latter Vandermonde determinant is in¬ 
variant, up to multiplication by sgn(cr), under the transformations a G &n of the 
n-element set {0,1 ,..., n — 1} of indices and that 

E ^Ul)^'(2) • • • 4(“n-l) = 

cr£6n 


is the same determinant ([5]). 


□ 
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3. Some applications of the criteria 


For a (real or complex) a 7^ 0,1, introduce the sequence of quantities 

= = n = 0.1.2 ..... (6) 

Applying the binomial theorem to each of the factors in the numerator of integrand 
and then integrating we hnd out that 


j,l=0 

n 

= E 

j,e=o 


n\ n 


(_l)n+i+£^n-£ / 


n\ n 


(- 1 ) 


n+j+l 




-— + (loga) V 
j+E-n 


^ / \ 2 
n 


a-’, 


where logo is understood as the integral of {dx)/x along a given path from 1 to a 
in (E]). Since \j + ^ — n\ < n m each summand of the hrst sum, we conclude from 
the representation obtained that 


e Z[a] logo + Z[a], 7 t, = 0,1,2,..., (7) 

where denotes the least common multiple of the integers from 1 to n (and do = !)• 

1 in 

Note that dn ^ e as n —)■ 00 by the prime number theorem. If a > 1 is an integer 
then the inclusions ([7]) simply mean that dnin £ ^ logo + Z for n = 0,1, 2,... . 
The family above corresponds to the choice 

, , {x — l){a — x) , , da; , r n ^ 

z(x) = -, a;(a;) = — and 7 = l,a C M, 

X X 

in the notation of (a)-(e). When a > 1, one easily hnds that 

max {z{x)} = z{^/a) = {y/a — 1 )^. 

l<x<a 


Theorem 1. log 2 and log 3 are irrational. 

Proof. As (V^ — l)^e = 0.4663... < 1, the irrationality of log 2 follows already 
from application of Proposition [T] In the case of log3 we use (\/3 — l)^e^/^/4 = 
0.6004 ... < 1 and Proposition [21 □ 


4. A NEW PROOF OF THE IRRATIONALITY OF TT 

We can also use the above argument for a = i = \/—l, when the integrals in ([6]) 
produce approximations to 7r/2 = —i log* with coefficients from Q[i]. (For all 
practical purposes we can think of integration in ([6]) as going along the arc of the 
unit circle.) In fact, the change of variable x = (1 + it)/(l — it) transforms the 
integrals into 

/•i — fY 

4 = 4(z) = 2”+^z(-l-z)”y^ n = 0,1,2,...; 
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therefore, it follows from ([7]) and the latter that 

2i(l — G Ztt + Z, n = 0,1,2,..., (8) 

in this case, where { • } denotes the fractional part. 

Theorem 2 . tt is irrational. 


Proof. We nse 


(x — l)(i — x) 
max - 

a;garc(l,i) X 


2^/2 max 
te[o,i] 


l + P 


2 - \/2 = 0.5857... 


and the inclnsions (jH]). Since (2 — \/2)e^/^/4 = 0.6563 ... < 1, (a slight adaptation 
of) Proposition [2] implies that the approximated nnmber vr is irrational. □ 


We remark that the potentials of the integral constrnction 




n+1 


dt e 


+ 


from the classical perspective (that is. Proposition [T]) are already discnssed in the sec¬ 
tion “A second attempt” in [21 p. 375]: The hrst few approximations to vr prodnced 
by the integrals look promising, “[njnfortnnately in the long run the asymptotics 
have decided otherwise.” Proposition [2] can be applied (though not directly since 
the corresponding there have to be replaced with y/2 for n odd to accommodate 
the arithmetic part of the argument) to the integrals by choosing 


z{t) = 


- 1 ) 


u{f) = 


dt 


7 =[0,1] 


l + f2 ’ 1 + ^2’ 

and 6n = ddn- Note the related elementary identity 

det (c^+^nj+£) = det {vj+i), c G C \ {0}, 


(9) 


that allows to easily manipulate with the extra factors like 2^/^ above in computing 
the Hankel-determinant asymptotics. 


5. Commentary 


Though Proposition [2] does not sound very practical, a question is about how 
much we can relax hypothesis (d) to still possess its implication. 

One possibility is to “density” the sequence of linear forms and consider instead 
the sequence = r\n/k\ for some hxed integral k >2 and together with the corre¬ 
sponding Hankel determinants Rn = deto<i,j<n(f^i+i)- Surprisingly enough (as not 
discovered in the existing literature) the underlying “fc-stuttering” Vandermonde 
determinants 


det 

o<i,r<n 


) 


(that replace the classical Vandermonde determinants ([5]) in the above proof of 
Heine’s identity (E])) factors into the product of powers of Zj and of k Vandermonde 
determinants; roughly speaking, Rn behaves similar to Rfn/k\y leading to no 
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refinement of Proposition [2l Furthermore, we recall that the non-vanishing of Rn 
for inhnitely many indices n can be shown without the reference to a generalization 
of Heine’s identity but with the use of the equality 

OD OO 

= (1 + Z + . . . -t- 

and of an old result of Kronecker (see [51 pp. 566-567] or [T^ Division 7, Prob¬ 
lem 24]): the Hankel determinants deto<y^<n-i(D-i-j) for n = 0,1,... eventually 
vanish if and only if the generating series '^nZ^ represents a rational function. 

Further variations are still possible, for example, considering other even powers of 
the Vandermonde determinants in ([2]) — this has a nice interpretation by means of 
the Hankel hyperdeterminants (see [9]) or replacing the Vandermonde determinants 
by other polynomials p{z) e Z[zi,...,Zn\ (and ensuring that the corresponding 
integrals over 7 ” do not vanish). But there is a natural analytical obstruction to 
getting anything weaker than eA^/^/4 < 1. It is based on the fact that the n-variate 
Z-trans£nite diameter tz([0, £]"■), which is introduced and studied in [3] for general 
sets E C C"", tends to e/i as n — >■ 00 — see Proposition [3] below. This means that 
for a nonzero n-variable polynomial p{z) of total degree N with integral coefficients 
we have 

sup \p{z)\>p^, 

z&[0,e]^ 

where —)■ e/4 as n —)■ 00 , thus showing that the upper estimate (| 1 ]) is best possible, 

as the total degree of the polynomial ni<j<£<n(% “ ~ 

Proposition 3. For the multivariate Z-transfinite diameter introduced in |3] we 
have 

lim tj^{[a,hY) = 
n^oo 4 

Proof. There are two C-extensions of the n-variate Z-trans£nite diameter, called 
t{E) and T{E) in [15]. The Hilbert-type relation between tz{E) and t{E) = tc{E) 
is established in [3l Theorem 3.1]: 

tz{E) < ( 10 ) 

The estimate 

T{E) < tz{E) (11) 

trivially follows from the dehnition of the two characteristics. Furthermore, the 
bound T{E) > c{E), where c{E) is the capacity of E is proved in [15]. If we 
restrict our attention to the cartesian product E = [a,b]'^ then c{E) = c{[a,b]) 
by m property fl')] and t{E) = r([a, 6 ]) by [TS] property fl)]. It remains to use 
r([a, 6 ]) = c([a, 6 ]) = |6 —a|/4 to conclude from flTOl) and flTTl) that the required limit 
relation holds true. □ 

Our further remark refers to the fact that the data 7 C z(x) and u(x) from 
condition (d) in Section [T] give rise to the scalar product on a space of single variable 
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functions of z G M: 

{u,v) = / u{z{x))v{z{x)) u{x), 

J-y 

so that the sequence in (d) is nothing but the sequence of the corresponding 
moments. It may be of some interest to study the orthogonal polynomials Pniz) 
for n = 0,1,... arising from the product: the polynomial Pn{z) can be explicitly 
written as the determinant Rn+i in ([3]), in which the last row is replaced with the 
row 1, z,..., z"", or as the integral 



Plenty of the theory of orthogonal polynomials naturally extends to these general 
settings, though it is not clear how much of this can be used in the irrationality 
context. 


6. Catalan and g-A pery constants 
It is worth mentioning that Proposition [2] is applicable to the integrals 


dz 


[ 0 , 1]2 


x{l — x)y{l — y)\^ dxdy 


xy 


e ZC(2) + z 


and 


dt 



[ 0 , 1]3 


x{l - x)y{l - y)z{l - z) 


xy 
dx dy dz 


e zc(3) + z 


1 — (1 — xy)z J 1 — (1 — xy)z 

used by Beukers in his proof [1] of Apery’s theorem about the irrationality of ^(2) 
and of C(3)) Ih® Apery constant. This follows from 


max 


x(l - x)y{l -y) [ Vd-l 


xy 


P,y,2)e[0,i]3 ^ 1 ™,. ~ 


1 — (1 — xy)z 
the prime number theorem and calculation 


4 V 2 


= 0.4527... < 1, 


-(\/2- = 0.6624. 


< 1 . 


For the related construction of rational approximations 


-)4n+lj2 
Rn 


to the Catalan constant 


[[ ( 

1 

1 

J J[0,1]^ ' 

V l-xy ) 


1 — xy 


G ZG + Z 




-1)* 


k'^o 

(see [nmaiin]), none of the irrationality criteria works, because of the too impetu¬ 
ous growth of the denominators required. 
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Other ways of “densifying” the sequence of rational approximations, different 
from the one discussed in Section El can be proposed, that take into account the 
expression of z{x). For example, in the case of the Catalan constant we can choose 
rational approximations to be 

^L("+1)/5J _ 2;)LO+2)/5j^LO+3)/5j^]_ _ ^^L("+4)/5j 

(1 



X 


X —y) ^/^dxdj/ 

1 — xy 


G Q G + Q, 


u = 0,1,2,..., 


because the corresponding Hankel determinants without the integer parts in the 
exponents behave asymptotically like (e^/^/4)"' , where £ = ((\/5 —1)/2)®. However, 
numerical computation of the honest Hankel determinants Rn = deto<y£<n(hj+£) 

suggests Rn"' (e/4)^/^ as n ^ oo] this numerical observation can be rigorously 
justihed using again the “5-stuttering” Vandermonde determinants. 

Another interesting question is about possible g-extensions of Proposition [21 the 
work [B] suggests that there may be some. For example, in [7] g-analogues of the 
Apery-Beukers rational approximations to C(3) were constructed, which approxi¬ 
mate the g-series 


Cq(3) = ^cT2(fc)g^ = 

k=l m=l 


TYi'^qm 

l-qm 


^ q^{l q^) 


for g the reciprocal of an integer different from 0 and ±1. It is expected that this 
g-Apery constant is irrational for such values of g, but the corresponding g-approx- 
imations from [7] do not produce any irrationality result, because |r,|V- ^ 1 
while ~ |g|-9/-' as n —)■ oo, where G Z^g(3) + Z. Though no arith¬ 

metic consequences come out from considering the Hankel determinants Rn = 
deto<j^i<n{rj+i), an analytic argument [T71 Section 4] shows their (better than ex¬ 
pected) behaviour |,|‘/3 as n —>■ cxo. 

More ideas are required to deal with the quantities like G and Cg(3)- 
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